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Abstract 

We consider dynamic network problems in which the control and disturbance inputs take values in 
finite sets. We derive a necessary and sufficient condition for the existence of robustly control invariant 
sets. We show that a stronger version of this condition is sufficient to guarantee robust global attractivity, 
and we construct a counterexample demonstrating that it is not necessary. Our proofs of sufficiency, being 
constructive, allow us to immediately extract the corresponding robust control laws as well as to establish 
the invariance of given sets. We conclude our study with two simple illustrative examples. 

1 Introduction 

Production networks, distribution networks and transportation networks are all instances of dynamic network 
problems in which the inventories, controlled flows and demands evolve in time. The dynamics of these 
systems are linear: The state of the system represents accumulated products in production stages, inventory 
stored in warehouses or commodities available at hubs. The control inputs represent the controlled flows 
between these production stages, warehouses or hubs, while the disturbance inputs represent fluctuating 
supplies of raw materials as well as unknown consumer demands. Practically, excess inventory results in 
undesirable holding costs, while shortages result in disruptions in the production or supply chain and a 
potential loss of clients. In view of this tradeoff, it is desirable to maintain inventory in all parts of the 
network within reasonable bounds in spite of the operating uncertainty. 

In this context, questions of stability, reachability and robustness naturally arise. Such questions have 
been previously considered in dynamic networks where the unknown but bounded disturbances as well as 
the control inputs are assumed to take values in continuous sets, and the invariant sets of interest range from 
hyperboxes to polytopes, ellipsoids, and more generally, convex sets [^ [TnHTB] . Related research directions in 
operations research address optimization problems in uncertain, multi-period inventory models. Tradition- 
ally, cost functions penalizing holding and storage costs are used in conjunction with dynamic programming 
techniques to robustly maintain the stock (i.e. the state) within a neighborhood of zero. Base stock policies 
(s, S), in which stock is reordered up to level S anytime it falls below a lower threshold s, can be interpreted 
in terms of guaranteeing upper bounds of robustly control invariant sets. Motivated by the desire to over- 
come the dimensionality problems inherent in dynamic programming and to mitigate the bullwhip effect and 
other undesirable dynamics, newer, robust optimization based approaches f8^ , including Affinely Adjustable 
Robust Counterpart (AARC) methods. Globalized Robust Counterpart (GRC) methods [5], and various 
extensions, have been developed. Often in the supply chain literature, the network topology is restricted to 
a tree, chain, or other special structure that can be exploited in the solution. 

Robustly control invariant sets have also been studied in a more general control context due to their 
relevance to a variety of control problems, including constrained control and robust model predictive control. 
Early results exploited dynamic programming techniques to characterize invariant sets [SJ[7]- An in-depth 
survey of the existing body of literature in control, emphasizing techniques that exploit connections with 
control Lyapunov functions, dynamic programing, and classical analytical results can be found in [3]. New 
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characterizations of robustly control invariant sets for discrete-time linear systems were recently proposed 
in |16| . allowing for tractable, convex-optimization based analysis approaches. 

We consider general discrete-time dynamic networks in which the inputs are restricted to take their values 
in finite alphabets, and we focus on a class of polytopic invariant sets (hyperboxes) . We seek to address 
the following two questions: Under what conditions does a robustly control invariant hyperbox exist, and 
under what conditions is such a set robustly globally attractive? Our contribution consists of two main 
results: First, we derive a necessary and sufficient condition for the existence of robustly control invariant 
hyperboxes. Second, we show that a stricter version of the above condition is sufficient, though not necessary 
in general, to guarantee robust global convergence of all state trajectories to the invariant set. While the 
problems are posed as existence questions, our results touch on analysis and synthesis problems as well. 
Indeed, the proposed proofs of sufficiency are constructive: As such they allow us to immediately extract 
the corresponding robust control laws. Moreover, they establish sufficient conditions for verifying invariance 
of a given set. 

The main novelty in our setup, which distinguishes it from the above referenced literature, is in the discrete 
nature of the control and disturbance inputs in conjunction with set based models of uncertainty. Discrete 
action spaces are justifiable from both practical and theoretical perspectives: From a practical standpoint, 
goods are usually processed, transported and distributed in batches. From a theoretical standpoint, the study 
of systems under discrete controls and disturbances has sparked much interest in recent years as evidenced 
by the literature on alphabet control [T31[T71[TH] , mixed integer model predictive control [T], discrete team 
theory [TH] and boolean control [3] where the problems of interest are often formulated as min-max games [5] . 
Moreover, it should be noted that in contrast to the supply-chain or multi-echelon, our work considers very 
generic network models. 

Organization: We present the problem statement in Section 12.11 and explain its practical significance in 
Section [2^ We state the main results in Section[3l present a complete proof in SectionlH and present simple 
illustrative examples in Section [5l 

Notation: M, Z, R+ and Z+ denote the reals, integers, non-negative reals and non-negative integers, 
respectively. [x\i denotes the i^^ component of a; S M". hull{S} and int{S) denote the convex hull and 
interior, respectively, of S* C R". B" denotes the set of vertices of the unit hypercube, that is B" = {0, 1}". 
For ^ C Z, M e Z"^™, MA"^ denotes the image of A"" by M, that is MA"" = {a £ Z"|a = Mb for some b e 
A""}. Given a set X = [0,x+] x . . . x [0, C M", Vx denotes its set of vertices, that is Vx ^ {x e M"|[a;]i e 

{0,4}}- 

2 Problem Setup 

2.1 Problem Statement 

Consider the dynamic network described by 

x{t + 1) = x{t) + Bu{t) - Dw{t) (1) 

where t G Z_|_, state x{t) E M", control input u{t) E and disturbance input w{t) g W^. The control and 
disturbance alphabet sets U — {ai, . . . , a^} C Z and W = . . . , bq} C Z, respectively, are ordered with 
ai < . . . < and 6i < . . . < feg. Be Z"^", D e Z"^p are given. 

Definition 1. A hyperbox X = [Q,x^] x . . . x [0,2:^] C M" is robustly control invariant if there exists a 
control law (p : X ^ U"^ such that for every x{t) E X , x{t + 1) — x{t) + Bip{x{t)) — Dw{t) E X for any 
disturbance w{t) E W^. 

Remark 1. When X = [0,x^] x ... x [0,a;^] is robustly control invariant, then so is any other hyperbox 
X' = [x^ , Xi + x^] X ... X [x^ , x^ + x^]. Indeed, control law ip' : X' — > defined by (p'{x) — (p{x — x^), 
where [x^]i — x^ , verifies this assertion. 
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Definition 2. A hyperbox X = [0, x^] x . . . x [0, x^] C R" is robustly globally attractive if there exists a control 
law ijj : M" \ X — >■ U"^ such that for every initial condition x{0) £ R" \ X and disturbance w : Z+ — > W^, the 
corresponding state trajectory satisfies x{t) G X for some r G Z+. 

We are interested in answering two questions for the dynamic networks described in ([T]): 
Question 1: Under what conditions does a robustly control invariant hyperbox exist? 
Question 2: Under what conditions does a robustly globally attractive and control invariant hyperbox 
exist? 



2.2 Relevance of the Model and Problem Statement 

The dynamics in ([T]) describe a unified but fairly general abstract model for three types of logistic networks: 
Production networks, distribution networks and transportation networks. In the first scenario, the nodes of 
the network represent "products", be they raw materials, intermediate products or finished products. The 
i*^ component of the state vector thus represents the amount of product i. The hyperarcs of the network 
represent production processes or activities, some of which may be fully or partially controlled by the 
operator of the network. Additionally, the network may interact with its external environment through both 
controlled and uncontrolled flows representing (generally uncertain) supply of raw material and demand of 
various products. The Bu — Dw term thus encodes the various production processes, supplies and demands, 
with matrices B and D representing the network topology and inputs u and w representing the controlled 
and uncontrolled flows, respectively. Likewise in the second and third scenarios, the nodes of the network 
represent warehouses and transportation hubs, respectively, with the i*^ component of the state vector thus 
representing the quantity of commodities present in the i*^ warehouse/hub. The Bu — Dw term encodes the 
various transportation routes, distribution protocols, supplies and demands, with matrices B, D representing 
the network topology and inputs u and w again respectively representing the controlled flows and uncertainty 
in the system. 

In this setting, it is intuitively desirable to contain each component of the state vector within two bounds, 
a zero lower bound and a positive upper bound. In the case of production networks, the lower bound guards 
against shortages and interruptions in the production process. In the case of distribution and transportation 
networks, the lower bound guards against the underuse of distribution and transportation resources. In 
all scenarios, the upper bound ensures that the storage capabilities of the system are not exceeded. The 
question of existence of robustly control invariant sets, specifically hyperboxes (Question 1 in Section l2.ip . 
thus naturally arises. 

Moreover in this setting, the model of uncertainty (specifically the choice of set W) encodes the typical 
uncertainty encountered in day to day operations. Since it is impossible to rule out rare occurences of large 
unmodeled uncertainty (be they emergencies or catastrophic events) that would drive the system away from 
its typical operation, it is reasonable to question whether the system can recover from such events: The 
question of robust global attractivity of the robustly control invariant hyperboxes (Question 2 in Section 
12. ip thus naturally arises. 



3 Main Results 

Consider the following sets for i G {1, . . . ,n}: 

U[ = {ueW\[Bu- Dw]^>oyw eWP}, 
Uf = {ueW\[Bu-Dw]^<OyweyVP}, 

Ul, = {u eZ^"|[Su-Z3w], < 0,Vw e WP}. 

Associate with every x G R" a signature, namely an n-tuple (si, . . . , s„) with = + if [x]i = and Sj 
if [x]i > 0, and two subsets ofU"^ defined by 

u,=ul^n...r\Ul^, 
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We can now state the main results. The first provides a complete answer to Question 1: 
Theorem 1. The following two statements are equivalent: 
(a) There exists a set X = [0,Xi] x . . . x [0,a;^] that is robustly control invariant, 
(h) The following condition holds 

U, 7^ 0, Vz e B". (2) 

The second proposes a sufficient condition for the existence of a robustly control invariant set that is 
globally attractive, thus giving a partial answer to Question 2. In general, this condition need not be 
necessary as discussed in detail in Section |4l 

Theorem 2. // the following condition holds 

U* 7^ 0, Vz e (3) 
there exists a robustly control invariant and globally attractive set X = [0,a;^] x . . . x [0,x+]. 

Remark 2. While it is somewhat diss appointing that both conditions ^ and ^ are combinatorial in 
nature, and thus grow exponentially with the number of nodes in the network, the result is not surprising. 
Indeed in \l(A\lS\j . discrete-time networks subject to analog control inputs and disturbances were analyzed. 
For instance, it was shown that when W = [d^,d'^], hi = for some scalars d^ < d'^ , < u'^ , the 

condition 

DWP C int{BW) 

is necessary and sufficient for the existence of a globally attractive robustly control invariant set. While 
this condition appears deceptively simple, verifying it is known to be NP-hard: Indeed, in general it requires 
checking 2" — 2 constraints \13}/ . 

4 Derivation of Main Results 

4.1 Existence of a Robustly Control Invariant Set 

We begin by establishing a necessary condition for a given control law to render a set robustly control 
invariant. 

Lemma 1. Let X = [0, x^] x . . . x [0, x^] be robustly control invariant, and consider a control law (p : X ^ U™ 
as in Definition[l[ Then at every vertex x G Vx, we have 

(p{x) e Uz 

where z is the unique element o/B" whose signature is identical to that of x. 

Proof. Assume that Lp{x) ^ lAz for some x e Vx with signature (si,...,S2), and consider a set X = 
[0,a;f] X ... X [0, a;+], x'l > 0. Pick awe \Uz. By assumption, there exists an i e {!,..., n} 
such that u ^ HI.. Now letting x{t) — x and applying control input (p{x){t) = u, we have [x{t + — 
[x{t) + Bu{t) - Dw{t)]i = [x]i + [Bu - Dw]i, and [x{t + - [x{t)]i = [Bu{t) - Dw{t)].i satisfies 

[x{t + 1)], - [x{t)\, > 

for some w G VV when [x\i and 

[x{t + l)], - [x{t)], < 

for some w G when [x]i = 0. Hence x{t + l) ^ X ioi some w G W^. Noting that the choice of m £ 

was arbitrary allows us to conclude that X is not robustly control invariant. □ 
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Remark 3. Lemma{l\can be be interpreted as the counterpart to the necessary condition in Nagumo's The- 
orem '151, adapted to the discrete-time, discrete alphabet setting of interest here. Nagumo's sub-tang entiality 
condition, while known to be sufficient for continuous-time systems and linear discrete-time systems under 
analog inputs, is not sufficient for general discrete-time systems. As such, it is not surprising that additional 
constraints need to be placed on the set to ensure sufficiency of condition (0) in establishing set invariance. 
This is presented next. 

Lemma 2. // the following condition holds 

U, 7^ 0, Vz e B" (2) 
the set X = [0, 2ii] x . . . x [0, 2i„] is robustly control invariant for any choice Uz £ Uz whenever Li > L° ~ 



max 



[Buz — Dw]i 



Proof. Assume that Uz ^ ^ holds for all z G B" and pick for each z G B" a control input Uz ^ Uz. The 
set X = [0, 2Li] X ... X [0, 2L„] is robustly control invariant. Indeed, consider the control law if : X ^ U"^ 
defined by 1^9(2;) = Uz{x), where z{x) G B" is the unique vertex of the unit hypercube with signature Si = + 
if [x]i < Li and — otherwise. Note that under this control law, we have for every i £ {1, . . . , n} 

[x{t + 1)], = [x{t)], + [B^{x{t)) - Dw{t)], 

Thus by construction, when < [x{t)].i < L„ < [Bip{x{t)) - Dwit)], < L° < L, and < [x{t + 1)], < 2Li. 
Likewise when L, < [x{t)]^ < 2L,, -L, < ~L° < [Bip{x{t)) - Dw{t)], < and < [x{t + 1)]^ < 2L,. It 
follows that X is robustly control invariant. □ 

Remark 4. Note that when the sufficient condition holds, our proof effectively provides an immediate 
construction of a full state feedback control law rendering any sufficiently large hyper box robustly control 
invariant. For any choice ofuz G Uz, z G B", the corresponding robustly control invariant set X = [0,2L°] x 
. . . X [0, 2L°] can be interpreted as an outer set (or superset) of the smallest robustly control invariant set. 

We are now ready to prove Theorem [T] 

Proof of Theorem\^ (a) (b): Assume that Uz = % for some z G B" with signature (si, . . . , s„), consider 
a set X = [0, x'^] X ... X [0, a;+], x^ > 0, and pick vertex x G 14 with signature (si, . . . , S2). For any control 
law (/J : X — >■ W™, we necessarily have '^{x) ^ Uz (as the latter set is empty). Noting that the choice of control 
law was arbitrary allows us to conclude that X is not robustly control invariant. Noting that the choice of 
X was arbitrary allows us to conclude that a robustly control invariant set cannot exist, thus completing 
our argument. 

(b) => (a): Follows directly from Lemma [2] which provides an explicit construction for such a set. I 
4.2 Robust Global Attractivity 

We begin by establishing sufficient conditions for a set to be robustly globally attractive. Our proof hinges 
on the construction on an appropriate Lyapunov-like function. 

Lemma 3. // the following condition holds 

U; 7^ 0, z G B" (3) 



set X = [0, 2LI] X ... X [0, 2L*^ with L* = max [Buz - Dw], 

Uz eU* 



is robustly globally attractive for any choice 



Proof. Pick a choice Uz eU*, for z G B". To prove global attractiveness of the corresponding set X, consider 
the control law ( : M" \ X ^ U™ defined by ((x) — Uz(x), where z{x) G B" is the unique vertex of the unit 

hypercube with signature = + if [x]i < Li and Si = — otherwise. Letting A = min [Buz — Dw]i , note 

that A > 0. Consider the function : M" ^ M defined by 

V{x) ~ max min [x]i — [y]i 
i yex 
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Note that by construction, we have V{x) > 0, for all x G M" and V{x) = x E X . Moreover, observe 
that under control law ^ we have 

V{x{t + l)) -V{x{t)) <0 (4) 

along system trajectories whenever x{t) e M" \ X. To verify this, consider any x{t) G M" \ X, and consider 
the i*'* coordinate direction. We have: 



mm 



< min 



(5) 



with equality holding only in the case where the right hand side is zero. Indeed, we can distinguish three 
cases: 

• [x{t)]i > 2L*: In this case, by construction we have < [x{t + 1)]^ < [a;(i)]i, and 



max{0, [x{t + 1)], - 2L*} = min [x{t + 1)], - [y], 

yex 



< mm 

yex 



[xm - [y]^ = - 2^: 



< < 2L*: In this case, by construction we have < [x{t + < 2L*, and 



mm 

vex 



[x{t + l)],-[y]. 



= mm 
yex 



= 



• [a;(t)]i < 0: In this case, by construction we have [2:(i)]i < [x{t + l)]i < 2L*, and 



max{0, -[x{t+ = min [x{t + ~ [y], 
v&x 



< mm 

vex 



Equation ([4]) thus follows from the the fact that ([5]) holds for each i G {1, . . . with equality only when 
both terms are identically 0. 
Moreover, we have 

V{x(t + l)) <V{x{t)) ~ A (6) 

whenever x{t) G R" \ X and x{t + 1) G M" \ X. 

Thus, for any choice of initial condition x(0) and of disturbance input w : Z-|_ W^, we conclude 
from (HI) that lim V{x{t)) 0. Moreover, we conclude from ([5]) that there must exist a r > such that 

V(x{t)) = 0, or equivalently x{t) E X. Hence X is globally attractive. □ 
We are now ready to prove Theorem [21 

Proof of Theorem\^ The proof is constructive. Assume that ([3]) holds and pick a choice E U*, for 
z G B". The set X — [0,2LJ] x ... x [0,2L*] is robustly control invariant by Lemma [21 since U* C Uz- 
Moreover, X is globally attractive by Lemma [31 ■ 

In the case of a degenerate network consisting of a single node (i.e. when n — 1, p and m arbitrary), 
condition ([3]) is necessary as well as sufficient for robust global attractivity: 



Proposition 1. When n = 1, if U* 
robustly globally attractive. 



for some z G {0, 1}, there cannot exist a set X — [0, L] that is 



Proof. Assume, without loss of generality, that Ul = % and consider a set AT = [0, L]. It follows that for 
any x{t) > L and any choice of control law, x{t + 1) > x{t) for some 'w{t) G W^, call this Wo- We can thus 
always find a disturbance input, namely w(t) = Wo-, i > 0, for which any state trajectory initialized to the 
right of the interval [0, L] will remain to the right of the interval for all times, and X is not robustly control 
invariant. □ 

In general, however, condition (jS]) is not necessary for global attractivity. Indeed, consider the second 
order counterexample constructed as follows: 
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Table 1: Payoff matrix for zero sum game in Example 1 



u/w 


-6 


4 


-100 


-100B + 6D 


-lOOB 


-AD 


-2 


-2B + 6D 


-2B- 


-AD 


3 


3B + 6D 


3B - 


AD 


150 


150B + 6D 


150B- 


-AD 



Counterexample 1. Let n = 2, m 

we have four possible control inputs, 

= |ui = 

Computing the relevant sets, we get 



2,p=l,W = {0}, U = {-1,3}, and B 



In this case. 



' -1 ■ 




■ -1 ■ 




CO 




■ 3 " 


-1 


,W2 = 


3 


,■"3 = 


-1 


,U4 = 


3 



U^ = {ui,U3}, Ul,-- 


= {wi}, 


U\ = {u2,U'i,Ui}, 




{U2,'"4}, 


Ul = {ui,M2,'"4}, Wi^ -- 


= {"2}, 


Ul = {mi,M3,U4}, 




{"3}. 


On the unit hypercube, we have 










^0,0]' 


= {U3,U4}, 


^[0,0]'=^' 






^0,1]' 


= {U2,U4}, 








^1,0]' 


= {U1,U3}, 








^14]' 


= {Ul}, 








By Lemma\^ there exists a set that is 


robustly control invariant, namely X 


= [0,24] 


x[0,8]. While 



(0) does not hold, it is easy to note that X is also robustly globally attractive. Indeed, consider the control 
law ip -.K? -^U"^ defined by 



Lp{x) 




when [x\i > 12,0 < [x\2 < 8 
when [x\i < 12,0 < [x\2 < 8 

when [x]2 < 

when [x]2 > 8 



It is straightforward to verify that (p renders X robustly control invariant and globally attractive. 



5 Illustrative Examples 

We begin with a simple scalar example for intuition, followed by an example of a production network 
consisting of six nodes and ten arcs. 

Example 1. Consider the scalar dynamics (n — m — p ~ 1) given by 

x{t + 1) = x{t) + Bu{t) - Dw{t) 

with alphabets U = { — 100, —2, 3, 150} and W = {—6, 4}. We begin by computing sets W-f and (no need 
for indices 'i' in this case) by inspecting Table\^ whose entries are simply the values of 'Bu — Dw': This 
table can be interpreted as the payoff matrix of a zero sum game between players u and w. 
We have and U- ^ IJ) iff 

(-100B + 6D<0 ( 150B + 6D>0 

\ -lOOB ~AD <0 \ 150B -AD>0 " 

We thus conclude that a robustly control invariant set indeed exists iff B > max{0.61?, — 0.04Z?} , and is 
moreover globally attractive provided strict equality holds. Consider for example the case where B = D = 1, 
for which a robustly control invariant set is guaranteed to exist: It is starighforward to verify that X — [0, 157] 
is robustly control invariant (the smallest such set!). 
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Example 2. We next analyze a production process described by the network with 6 nodes and 10 hyperarcs 
depicted in Fig. [7J The process produces four different products, AAB, ABB, AAAB, and ABBB (nodes 3 
to 6) from raw materials A and B (nodes 1, 2). Disturbances wi and W2 model the uncertain supply in raw 
material, while disturbances to wq model the uncertain demands for the various products. Control inputs 
ui and U2 model the nominal supply of raw materials A and B, while U3 to uiq model the various production 
schemes that can be utilized: For instance, product ABB can be produced by combining raw materials A and 
B (hyperarc U4 connecting nodes 1 and 2 to node 4), or from product AAB by releasing one unit of A and 
adding one unit of B (hyperarc uj connecting nodes 2-3 to nodes 1-4)- The quantitative relations among 
products and materials are captured by matrix 

1 -3 -1 " 
1 -1 -3 
1 
10 0' 
1 
1 

The perturbation in supply of raw material and the unknown demand enter the model as an additive uncer- 
tainty, with matrix D being the 6-dimensional identity matrix. We assume that U = {a Cz Z|0 < a < 400} 
and W = {6 G Z|20 < b < 40}. Under these assumptions, condition (0) holds. We proceed to verify this 
without explicitly constructing the sets U* , by employing a heuristic approach. Indeed, we solve the following 
linear program with decision variables X and u for each vertex z G ; 

min A 

subject to [Bu]i - 40 > A if = 
[Bu]^ - 20 < -A if [z], = 1 
[u], > 
[u]i < 400 
A > e 

where e > is a chosen parameter. Note that the linear program returns a fractional solution that can be 
rounded to an integer solution within the admissible input set. Also note that the existence, for each choice 
of z E B^, of an integer solution satisfying the first four LP constraints (for any A > Oj effectively ensures 
satisfaction of since these constraints represent extremal (worst case) values of the disturbance inputs. 
Parameter e allows us some flexibility in applying this heursitic: A larger e translates into a higher likelihood 
that the rounded integer solution will satisfy the desired constraints, at the expense of missing potential 
solutions if e is too large. Once we identify rounded integer controls and double check that condition ^ still 
holds, we store these values in a look-up table ready to be implemented in feedback form, chosen in agreement 
with the control law proposed in the proof of Lemma O Having chosen the feedback law, we can now also 
compute the robustly globally attractive control invariant set [0,2ij*] x ... x [0, 2Lg]; The obtained values 
are displayed in Table\^ 

Having computed the feedback control law offline, we report on Monte Carlo simulations of the closed loop 
system. We ran 30 different paths with 600 samples (horizon steps) each starting from uniformly randomly 



B = 



1 





-2 


-1 


-1 





-1 





1 


-1 


-2 





-1 


1 








1 





-1 





1 











1 





-1 


-1 














1 























1 
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Table 2: Size of the invariant box [0, 2L|] x . . . x [0, 2Lg]. 



LI = 625 


LI = 625 


LI = 221 


LI = 221 


Lt, = 237 L* = 237 



1000 
500 




1500 
1000 
500 



O P--'-. . 

200 400 600 200 400 600 



time 



time 



1000 
500 





200 400 600 200 400 600 



time 



time 



400 
200 




l]l Li,.i.ii.l.,Li-iilllnjJL^.a.,J i].t. 



1000 
500 



200 400 600 200 400 600 

time time 



Figure 2: First sample path of the Monte Carlo simulations for Example [2j 



selected initial states in the interval [— 1000,2L* + 1000] and with random demand uniformly drawn from 
the admissible intervals. The first sample path starting from initial state (1242, 1543, 1281, 779, —161, 1468) 
is depicted in Fig. [H with the dashed lines delineating the invariant hyperbox. All simulations are carried 
out with MATLAB on an Intel(R) Core(TM)2 Duo CPU P84OO at 2.27 GHz and 3GB of RAM. The run 
time of the offline computation of the control law is less than 5 seconds, while the run time of loading all the 
controls in the look-up table and running the Monte Carlo simulations is about 12 seconds. 



6 Conclusions &; Future Work 

We considered dynamic networks in which the control and disturbance inputs take their values in finite sets. 
We established a necessary and sufhcient condition for the existence of robustly control invariant hyperboxes, 
and we showed that a stronger version of this condition is sufficient, albeit not necessary, to guarantee robust 
global attractivity. The proposed conditions are combinatorial in nature, not suprising as the problem is 
known to be NP-hard even when the input signals are analog. 

Future work will focus on getting a handle on the size of the smallest such invariant hyperboxes, as well 
as considering more interesting models of finite alphabet uncertainty. 
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